Application of Nonlinear Conductivity Spectroscopy to Ion Transport in Solid Electrolytes 
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The field-dependent ion transport in thin samples of different glasses is characterised by means of nonlinear 
conductivity spectroscopy. AC electric fields with strengths up to 77 kV/cm are applied to the samples, and 
the Fourier components of the current spectra are analysed. In the dc conductivity regime and in the transition 
region to the dispersive conductivity, higher harmonics in the current spectra are detected, which provide infor- 
mation about higher-order conductivity coefficients. Our method ensures that these higher-order conductivity 
coefficients are exclusively governed by field-dependent ion transport and are not influenced by Joule heating 
effects. We use the low-field dc conductivity dijc snd the higher-order dc conductivity coefficient O-i^dc to 
calculate apparent jump distances for the mobile ions, aapp. Over a temperature range from 283 K to 353 K, we 
obtain values for agpp between 39 A and 55 A. For all glasses, we find a weak decrease of a^pp with increasing 
temperature. Remarkably, the apparent jump distances calculated from our data are considerably larger than 
typical values published in the literature for various ion conducting glasses. These values were obtained by 
applying dc electric fields. Our results provide clear evidence that the equation used in the literature to calculate 
the apparent jump distances does not provide an adequate physical description of field-dependent ion transport. 
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1. Introduction 

An important prerequisite for the application of solid ionic 
conductors as electrolytes in batteries, fuel cells and elec- 
trochromic devices is the preparation of materials with low 
electrical resistance. This can be achieved by either improv- 
ing the specific ionic conductivities of materials or by prepar- 
ing thin films. A rule of thumb for the technical applica- 
bility of solid electrolytes is an ionic conductivity exceeding 
10^^ S/cm. However, in the case of thin-film materials with 
thicknesses of 1 /im and less, conductivities of the order of 
10^^ S/cm to 10^^ S/cm can be sufficient. Thus, by means 
of thin-film technology, the range of solid electrolyte materi- 
als being suitable for technical applications can be extended 
considerably. In the case of crystalline materials, epitaxial 
methods are often used for the preparation of thin films, while 
sputtering methods and spin coating techniques are suitable 
for preparing thin films of glassy and polymeric electrolytes. 

When the ionic conductivity of a solid electrolyte layer does 
not depend on its thickness d, the electrical resistance of the 
layer is simply proportional to d. However, in the case of thin 
samples, the assumption of a thickness-independent conduc- 
tivity is often not valid. The reason for this is the presence of 
large electric fields in the samples, which lead to nonlinear ion 
transport effects. Let us consider an example: A film with a 
thickness of d= 100 nm is used in a lithium battery providing 
a voltage of f/ = 5 V. In this case, the electric field strength in 
the film isE = U /d = 500 kV/cm. It is well known that at field 
strengths £ > 50 kV/cm, the ionic conductivity of many solid 
electrolytes increases with increasing field strength. This in- 
crease of the ionic conductivity should generally improve the 
applicability of thin films. 

A possible mathematical description of the nonlinear elec- 
trical properties of ionic conductors is the following power se- 
ries for the dependence of the current density j on the electric 
field strength E: 



are higher-order conductivity coefficients. Eq. ([0 contains 
exclusively odd terms, since the function j{E) is an odd func- 
tion. 

A number of field-dependent conductivity studies have 
been caiTied out on ion conducting glasses (lll2.l3.BS01- In 
these studies, dc electric fields, Edc, were applied, and the dc 
currents densities, j^c, were determined. Often the jdc{Edc) 
curves can be reasonably well fitted by a hyberbolic sine func- 
tion: 
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j = (7i£' + C73£'^ + C75-£'^ 
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Here, Gi denotes the low-field conductivity, while a^, G5 etc. 



Here, q denotes the charge of the mobile ions, while kg and 
T are Boltzmann's constant and the temperature, respectively. 
The quantity Oapp has the unit of a distance, and is, in the fol- 
lowing, called an 'apparent jump distance'. Eq. (|2ji can be de- 
rived theoretically in the framework of a random walk theory 
with mobile ions carrying out thermally activated hops in a pe- 
riodic potential landscape. In this framework, aapp is identical 
to the actual jump distance of the mobile ions, i.e. to the dis- 
tance between neighboring sites. However, the apparent jump 
distances derived from fits of experimental data using Eq. (|2|l 
are generally between 15 A and 30 AQl HH H H |i] . Thus, 
flapp is much larger than typical distances between neighbor- 
ing ionic sites in glasses. In the framework of molecular dy- 
namics simulations, these typical distances were found to be 
about 2.5 - 3 A|7, 8]. The physical origin of the large values 
for flapp is most likely related to the amorphous structure of 
the matrix the ions are moving in. This amorphous structure 
should result in a highly disorderd potential landscape with 
spatially varying site energies and barriers, in contrast to the 
periodic potential landscape used in the random walk theory 
mentioned above. However, up to now, there is no comprehen- 
sive theory relating the the shape of the potential landscape to 
nonlinear ion transport. 

A drawback of nonlinear conductivity studies by means of 
dc electric fields is that no direct experimental information 
can be obtained on Joule heating effects in the samples. Joule 
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heating leads to an increase of the sample temperature, which 
leads again to an increase of the ionic conductivity. This con- 
ductivity enhancement may pretend a nonlinear ion transport 
effect. We would like to emphasize that the differentiation be- 
tween nonlinear ion transport and Joule heating is very impor- 
tant for technical applications. While an increase of the ionic 
conductivity due to nonlinear transport effects should gener- 
ally improve the applicability of materials. Joule heating often 
leads to an electrical breakdown of materials. 

A differentiation between nonlinear ion transport and Joule 
heating is possible when ac electric fields are used. According 
to Eq. Q, the application of a sinusoidal electric field E{t) = 
Eq ■ sin(a)f) leads to the following expression for the current 
density being in phase with the electric field, / : 

/ = (7i'-£'o-sin(a)f) + (73-(£'o)^-sin^(cof) (3) 
+(75-(£'o)^-sin^(C(Jf) + 

= c[{(o) -fio ■sin(a)f) + - (73(0)) • (£0)^ •sin(cof) 

-i(7;^(3a))-(£o)^-sin(3a)r) 

+l^a^(c.).(/<o)^.sin(c.O 

--^ (7^(3«)-(£o)^-sin(3«f) 

+ -^- (y'^{5o}) ■ (£0)^ •sin(5«f ) + .... 
16 

An analoguous expression can be written down for the out- 
of -phase current density /' by replacing c7/((b) by (7/'(a)) and 
by replacing the sine functions by cosine functions. However, 
in the following, we show exclusively results for the in-phase 
current density /. As seen from Eq. (|3}, the third-order term 
in Eq. Q generates a Fourier component in the current den- 
sity spectrum at 3 ©. Accordingly, the fifth-order term gener- 
ates a 5 0) component, etc. Thus, the higher harmonic contri- 
butions can be used to determine the higher-order conductivity 
coefficients 03, 05 etc. In contrast. Joule heating of the sam- 
ple leads to an increase of the low-field conductivity a\, but 
does not generate higher harmonics in the current spectrum. 

Tajitsu used this ac method to determine the higher-order 
conductivity coefficients of some polymer electrolytes based 
on polyethlyenoxide / lithium salt mixtures |9, 10]. These 
measurements were done in a frequency range from 10-2 Hz 
to 10^ Hz. At low frequencies, the low-field conductivity as 
well as the higher-order conductivity coefficients become in- 
dependent of frequency and identical to their dc values. By 
using Eq. Tajitsu calculated apparent jump distances Oapp 
from these dc values, and he obtained values for flapp between 
40 A and 50 A. A drawback of Tajitsu's measurements, how- 
ever, is the usage of ion blocking metal electrodes. In this 
case, it is difficult to assess the respective influences on the 
higher-order conductivity spectra exerted by electrochemical 
processes at the interfaces between sample and electrodes, and 
by electron injection into the sample. 

Therefore, our approach is the utilisation of ion conducting 
liquid electrodes in combination with ac electric fields and a 
Fourier analysis of the current density spectra. As liquid elec- 



trodes, we use highly conducting salt solutions containing the 
same type of cations, which are mobile in the solid electrolyte. 
In this case, an electric field simply drives these cations from 
one liquid electrode through the solid electrolyte into the other 
liquid electrode. 

By means of this method we carry out systematic 
temperature- and composition-dependent studies on nonlin- 
ear ion transport effects in solid electrolytes. Very little is 
known about the respective influences on nonlinear ion trans- 
port exerted by the number density of mobile ions and by the 
structure of the solid matrix. In this paper, we present results 
on the nonlinear conductivity of three different ion conducting 
glasses. The compositions of these glasses are 0.127 Na20 
■ 0.096 CaO • 0.062 AI2O3 • 0.715 Si02 (NCAS12), 0.25 
NajO ■ 0.096 CaO • 0.062 AI2O3 • 0.592 Si02 (NCAS25), 
and 0.2 Na20 ■ 0.8 Si02 (NS20). The composition of the first 
glass is similar to Thuringian glass 0.101 Na20 • 0.084 CaO 
• 0.03 AI2O3 • 0.785 Si02 (+ trace amounts of K2O, MgO 
and Fe203), but the sodium oxide content is slightly higher. 
We chose this slightly higher sodium oxide content, since 
the room temperature dc conductivity of Thuringian glass is 
too low to be measured within our temperature and frequency 
window. For Thuringian glass andfora0.18Na2O • 0.82 Si02 
glass with a composition similar to our NS20 glass, values for 
the apparent jump distance flapp obtained from dc electric field 
measurements are available in the literature 1 1, 2]. 

By means of high-precision cutting and grinding tech- 
niques, samples with thickness between 65 jUm and 85/im 
were prepared. Much thinner samples are not suitable for 
our measurements, since the films have to be freestanding be- 
tween the liquid electrodes. However, since we are able to 
apply ac voltages with amplitudes up to 500 V, the resulting 
maximum electric field amplitudes of about 77 kV/cm are suf- 
ficient to detect higher harmonics at 3 a) in the current den- 
sity spectrum. Our measurements were carried out over a fre- 
quency range from 10 mHz to 10 kHz and over a temperature 
range from 283 K to 353 K. 



2. Experimental 

2.1. Preparation of thin glass samples 

Bulk glass samples were prepared from reagent grade 
chemicals of Na2C03, CaC03, AI2O3 and Si02. All chem- 
icals were mixed thoroughly and then melted at 1873 K for 5 
hours. After complete homogenization, the melts were poured 
into preheated stainless steel molds with a cylindrical shape. 
The obtained bulk glass samples were then annealed 40 K be- 
low their respective glass transition temperatures for 5 hours. 

Subsequently, the bulk samples were cut into cylindrical 
slices with a thickness of about 1 mm using a high-precision 
cutting machine (Struers Accutom-5) equipped with a dia- 
mond saw blade. The thickness of the slices was then further 
reduced by high-precision grinding using a lapping machine 
(Logitec PM5). Thereby, we attained sample thicknesses in 
the range from 65 /xm to 85 pLm. The maximum error in the 
thickness over the faces of a sample was about 2 jUm. 

The thin samples were attached to a highly resistive quartz 
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glass tube (conductivity < 10^'^ S/cm) by means of a high- 
voltage resistant Araldite glue (Vantico). The quartz tube with 
the attached sample was placed inside a quartz glass container, 
see Fig. [2 




■ - . , Highly resi stive 
quartz glass 

Araldite glue 
(high-voltage resistant) 



Thin glass sample 
Aqueous NaCl solution (1 mol/1) 



FIG. 1 : Schematic illustration of the experimental setup for nonlinear 
conductivity spectroscopy on thin glass samples 



Both the quartz glass tube and the quartz glass container were 
filled with a 1 molar aqueous NaCl solution, prepared with 
destilled water. Platinum wires connected to the high-voltage 
measurement system were dipped into the NaCl solutions. 
Since the electrical resistance of the NaCl solution is many or- 
ders of magnitude lower than the resistance of the glass sam- 
ples, there is virtually no voltage drop in the NaCl solutions, 
but the voltage applied to the platinum wires drops completely 
at the glass samples. Therefore, the NaCl solutions act as non- 
blocking electrodes. 

2.2. Nonlinear conductivity measurements 

Our high-voltage measurement system is based on the 
Novocontrol a-S high resolution dielectric analyser, which 
is equipped with a broadband high-voltage amplifier, a broad- 
band dielectric converter and a high-voltage sample head. A 
schematic illustration of the setup is shown in Fig. |2] The 
high-voltage system provides a frequency range from 3 jxHz 
to 10 kHz, a maximum voltage of 500 V, and a cuiTent resolu- 
tion of 5 fA. 

The digital waveforms of sample voltage and current are 
numerically Fourier transformed by a digital processor in the 
a-S analyser. Thereby, the amplitudes and phases of the base 
wave and of higher harmonics in the cuiTent spectrum are de- 
termined with respect to the sinusoidal voltage. Higher har- 
monics are detectable, if their amplitude is at least 10^^ of the 
base wave amplitude. 

The sample temperature was controled by the Novocontrol 
Quatro Cryosystem. 




FIG. 2: Schematic illustration of the Novocontrol high- voltage mea- 
surement system. 

3. Results 

In Fig.|3lwe show results for the base current density /(v) 
(closed symbols) and for the higher harmonic current density 
/(3v) (open symbols), measured after applying a sinusoidal 
electric field with amplitude Eq = 58.8 kV/cm to the NS20 
glass. Both Fourier components of the current density spec- 
trum are normalised by the field amplitude Eq and are plotted 
versus frequency v at three different temperatures, T = 293 K, 
313 K, and 333 K. The values obtained for /(3v) are negative. 

1 ^ 
According to the term — - • a^{3(o) ■ (Eq) ■ sin(3a)f) in Eq. 

a negative value for /(3v) indicates a positive value for 
(73 (3v). In order to include the /(3v)/£'o data in the log-log 
plot of Fig. |3l we have multiplied these data by a factor of -4. 
The f{v)/Eo curves reflect the frequency dependence of the 
low-field conductivity <7i'(v). At low frequencies, the low- 
field conductivity is independent of frequency and identical to 
the low-field dc conductivity. However, at higher frequencies, 
there is a transition into a dispersive regime with the conduc- 
tivity increasing with increasing frequency. With increasing 
temperature, the low-field dc conductivity increases, and the 
transition region from the dc conductivity to the dispersive 
conductivity shifts to higher frequencies. These frequency- 
and temperature-dependent features are well known for many 
solid ionic conductors. 

In the dc conductivity regime and in the transition region to 
the dispersive conductivity, the higher harmonic component 
j{3v) was detectable, see Fig.|3l However, at higher frequen- 
cies in the dispersive regime, the ratio /(3v)/ /(v) became 
so small, that reliable values for /(3 v) could not be obtained. 
The field dependence of /(3 v) at low frequencies was further 
analysed in order to calculate the higher-order conductivity 
coefficient CT^{v). According to Eq. (|3j, we can write the 
following expression for ./'(3 v): 



(4) 
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FIG. 3: Base current density /(v) (closed symbols) and higher har- 
monic current density /(3v) (open symbols), measured after apply- 
ing a sinusoidal electric field with amplitude Eq = 58.8 kV/cm to the 
NS20 glass. Both Fourier components of the current density are nor- 
malised by the field amplitude £() and are plotted versus frequency 
V at three different temperatures, T = 293 K, 313 K, and 333 K. In 
addition, the quantity /(3v) /Eq is multiplied by a factor of -4. 



Since at the applied field strengths, Fourier components at 5(0 
in the current spectrum were not detectable, the second term in 
Eq. must be considerably smaller than the first one. Thus, 
when we plot —4- /(3v)/£o versus (^o)^- the data should, 
in a good approximation, lie on a straight line with a slope of 
03 (3v). As an example we show in Fig. I^such a plot for the 
NS20 and NCAS25 glasses at a temperature of T = 293 K. 
For this plot, the frequencies were chosen such that /(v) and 
/(3v) are identical to their dc values. Note that in the case of 
the abscissa, we have normalised the amplitude of the electric 
field by the value 10^ V/cm. 

A closer look at Fig. [^reveals that the data do not exactly 
lie on a straight line, but on a curve with a slightly negative 

curvature. From this we conclude that the term ^ ■ O5 • (Eo)'^ 
in Eq. (|4} is not negligible. The solid lines denote fits of the 
data by a second-order polynomial. From the linear term in 
this polynomial, we obtain directly <7s dc- The quadratic term 
is negative, implying a negative value of 05 ij^-. Such negative 
values of (Jg ^/f we find for all glasses and at all temperatures. 
However, the error in the quadratic terms is so large that we 
can exclusively determine the sign of Os^^c, but not accurate 
values. 



4. Discussion 

Our results show that the field strengths applied in our ex- 
periments are sufficient to obtain accurate values for 03 ,./^. 



FIG. 4: Plot of -4 • /(3v) /Eq versus {Ea f for the glasses NS20 and 
NCAS25 at a temperature ofT = 293 K. The measurement errors are 
indicated by the size of the symbols. The solid lines denote fits of the 
data by a second-order polynomial. The linear and the second-order 
term of this polynomial provide information about O^jc and O^jc^ 
respectively. 

Additionally, we are able to determine the sign of C75 ^c- In 
the following, we first compare our results to literature results 
obtained by applying dc electric fields. As mentioned in the 
introduction, the jdc{Edc) data of many ion conducting glasses 
can, in a first approximation, be fitted by means of Eq. (|2j- A 
Taylor expansion of this equation in terms of the electric field 
yields the following expression for the ratio Oi^dcl (^i^c'- 

c>\,dc 6 ylksT J 

Solving this equation for the apparent jump distance aappj we 
arrive at: 

1 6 CJ3,/, [IkBTf 

flapp = -1 / o (O) 

This equation was used to calculate Oapp for the glasses. Our 
measurement technique ensures that these values reflect ex- 
clusively nonlinear ion transport effects and are not influenced 
by Joule heating. In Fig. |5j Oapp is plotted versus temperature 
T for the three glasses. Note that both axes are logarithmic 
axes. As seen from the figure, the two glasses NCAS12 and 
NCAS25 exhibit similar values of aapp- At room temperature, 
the apparent jump distance is about 45 A and decreases to 
about 40 A at temperatures in the range from 343 K to 353 K. 
However, in the case of the NCAS12 glass, the temperature 
dependence of a^pp is slighthly stronger. The data in this log- 
log plot can, in a first approximation, be fitted, by a straight 
line with a negative slope —a. This corresponds to a power 
law dependence of the apparent jump distance on temperature: 



5 



aapp °^ T~". For the NCAS12 glass we find a = 0.75 ± 0.18 
as compared to a = 0.51 ± 0.13 in the case of the NCAS25 
glass. The nonlinear ion transport in the NS20 glass is charac- 
terised by larger apparent jump distances. Oapp decreases from 
about 55 A at room temperature to about 50 A at T = 343 K. 
The temperature dependence is given by a = 0.5 ± 0.15 and 
is thus similar to that of the NCAS25 glass. 

60 ti I I I I I I I I 

: □ NS20: slope = -0.50 ±0.15 : 

: o NCAS12: slope = -0.75 ±0.18 : 




280 290 300 310 320 330 340 350 360 
T/K 

FIG. 5: Plot of the apparent jump distance Oapp versus temperature 
T for three different glasses (log-log plot). The solid lines denote fits 
of the data to the function aapp = T^". 



Remarkably, the apparent jump distances Oapp we have ob- 
tained from the ratio of (J^^dc/^^i.dc ^re even larger than typ- 
ical values published in the literature, which were obtained 
by fitting idc{Edc) data to Eq. (|2}. For instance, in the case 
of Thuringian glass with a composition close to our NCAS12 
glass, the literature value for the apparent jump distance is 
15.5 A ZagarandPapanikolauobtainedavalueof28.3 A 
for a 0. 18 Na20 • 0.82 Si02 glass with a composition close to 
our NS20 glass 12). The reason for this apparent discrepancy 
between our results and literature results is related to the va- 
lidity of Eq. (|2}. The fact that we find negative values for 
O^ jc for all glasses shows that Eq. (|3 does not provide an 
adequate physical description of field-dependent ion transport 
in our glasses. The validity of the hyperbolic sine function 
would imply positive values of all higher-order conductivity 
coefficients, including O^ jc- Since this is not the case, it is 
plausible that the determination of fl^pp from Eqs. (|3 and (|6} 
results in different values for Oapp. 

The electric field strengths applied in our study were not 
sufficient to detect harmonics in the current spectra at 5(0. 
These harmonics could be used to obtain accurate values for 
Os.i/f. As already mentioned, the second-order polynomial fits 
of the data in Fig. |4]yield only estimates for Os.dc- However, 
from these estimated values, we can pre-estimate that a reduc- 
tion of the sample thickness by a factor of two should lead to 
detectable harmonics at 5a). At the same time, this reduction 
of the sample thickness should allow a further extension of the 



frequency range where harmonics at 3 O can be detected. Such 
a reduction of the thickness should be achievable by means of 
chemical etching techniques. 

From a theoretical point of view, it is a great challenge to 
obtain a better understanding of the large apparent jump dis- 
tances aapp, of their temperature dependence, and of the neg- 
ative sign of O^ dc- We would like to emphasize that an ap- 
parent jump distance of the order of 50 A does not imply that 
ions perform single jumps over such a large distance. In the 
dc conductivity regime, we do not probe local hops between 
neighboring sites, but we probe the long-range ion transport. 
In the case of a disordered potential landscape with site en- 
ergy disorder, the distances relevant for the long-range ion 
transport could be distances between low-energy sites in the 
landscape and/or distances between sites close to the percola- 
tion energy. Detailed theoretical studies on particles moving 
in disordered potential landscapes will be required in order to 
obtain a detailed microscopic understanding of these relevant 
distances and their influence on the apparent jump distance 
flapp- Such studies are currently underway fTlll . 



5. Conclusions 

We have measured the nonlinear ionic conductivity of three 
different glasses in a temperature range from 283 K to 353 K. 
In the dc conductivity regime and in the transition region 
to the dispersive conductivity, the application of ac electric 
fields with amplitudes up to 77 kV/cm leads to harmonics 
in the current spectrum at 3 a), which can be used to calcu- 
late the higher-order conductivity coefficient a'^. From the 
ratio (J^^dc/ <^i.dc^ we derive values for the apparent jump dis- 
tance flapp of the mobile ions. In the temperature range of 
our study, we find values between 39 A and 55 A. The 
temperature dependence of the apparent jump distance is, in 
a first approximation, given by flapp °= T^", with a values 
ranging from 0.5 to 0.75. Our values for flapp are signif- 
icantly higher than literature values for various glasses ob- 
tained by fitting jdc{Edc) data to the hyperbolic sine func- 
tion jdc °= sinh(^flapp£'rft/(2feBr)). One reason for this ap- 
parent discrepancy is the fact that the hyperbolic sine function 
does not provide an adequate physical description of field- 
dependent ion transport. This can directly be concluded from 
the negative sign of Gs dc, which we find for all glasses. We 
pre-estimate that a further reduction of the sample thickness 
by a factor of two, which should be achievable by means of 
chemical etching techniques, will allow us to considerably ex- 
tent the frequency range, where harmonics at 3 a) can be de- 
tected, and, in addition, to obtain accurate values for o^^dc- 
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